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We report new limits on ultralight scalar dark matter (DM) with dilaton-like couplings to photons 
that can induce oscillations in the fine-structure constant a. Atomic dysprosium exhibits an elec¬ 
tronic structure with two nearly degenerate levels whose energy splitting is sensitive to changes in 
a. Spectroscopy data for two isotopes of dysprosium over a two-year span is analyzed for coherent 
oscillations with angular frequencies below 1 rad s _1 . No signal consistent with a DM coupling is 
identified, leading to new constraints on dilaton-like photon couplings over a wide mass range. Un¬ 
der the assumption that the scalar field comprises all of the DM, our limits on the coupling exceed 
those from equivalence-principle tests by up to 4 orders of magnitude for masses below 3 • 10“ 18 eV. 
Excess oscillatory power, inconsistent with fine-structure variation, is detected in a control channel, 
and is likely due to a systematic effect. Our atomic spectroscopy limits on DM are the first of their 
kind, and leave substantial room for improvement with state-of-the-art atomic clocks. 


Dark matter (DM) makes up the majority of matter 
density in our Universe. Its ubiquitous abundance can 
be measured through its gravitational influence, but lit¬ 
tle is known about the microphysical properties of the 
DM particle(s), such as the mass, spin, and any nongrav- 
itational interactions. If the DM is bosonic rather than 
fermionic, it can have a sub-eV mass and such high occu¬ 
pation numbers that it acts more like a classical wave— 
with frequency equal to its mass—than a particle. Light 
bosonic dark matter has a natural production mecha¬ 
nism, namely early-universe misalignment of the field rel¬ 
ative to the minimum of its potential w- Several mo¬ 
tivated candidates in this category exist in the literature, 
most notably the QCD axion [Mj and other axion-like 
particles, which, as parity-odd pseudo-Nambu-Goldstone 
bosons (PNGBs) of compact symmetry groups, primar¬ 
ily have derivative interactions with matter [7 . Light, 
parity-even bosons may also arise as PNGBs of noncom¬ 
pact groups such as those of scale, conformal, or shift 
symmetries, the most famous examples of which are dila- 
tons [SHE!. Small explicit breakings of these symme¬ 
tries may induce nonderivative operators for the scalar 
fields, such as mass terms and higher-dimensional oper¬ 
ators coupling them to matter. 

We focus on ultralight scalar fields <p with couplings to 
the (square of the) electromagnetic field tensor F^ v \ 

£ =5 \rnl4> 2 + 1 (1) 

where k = \Z4ttGn, Gn is Newton’s constant, and 
e ~ 0.303 is the electromagnetic gauge coupling (we use 
units in which H = c = 1). The interaction is normalized 
such that d e = 1 yields an attractive force of gravita¬ 
tional strength between electromagnetic energy densities 
at distances smaller than the inverse mass (r < mj) 1 ) 
through scalar <p exchange m Couplings d e 1 gener- 
ically arise if quantum gravity effects weakly break an 


underlying global symmetry of cj> near the Planck scale 
~ Gjy 1 ^ 2 [T2J. Quantum corrections proportional to d 2 
and naturalness considerations together suggest a min¬ 
imum mass-squared for cj>. However, given the existing 
hierarchy problems of the Standard Model, we remain 
agnostic to this issue and consider the full m^-de param¬ 
eter space (see Ref. m for more discussion). 

Equivalence principle (EP) tests such as the Eot-Wash 
experiment Il4j and Lunar Laser Ranging m constrain 
the coupling d e to be much less than unity: \d e \ < 3.6 • 
10~ 4 at 95% confidence level (CL) for < 3-10~ 14 eV. 
This limit will likely be improved to \d e \ < 7.8• 10 -5 with 
atom-interferometry techniques m ■ The EP-violating 
force that 4> mediates scales as \d e \ 2 , making vast im¬ 
provements to these limits challenging. 

Light scalar fields do not behave as perfect cold DM 
on short length scales, where their density perturbations 
have a nonzero sound speed. For < 10 -22 eV, they 
would have inhibited cosmological structure growth m- 
fl9l in conflict with observations, though these bounds 
disappear if <f> only makes up a small fraction of the dark 
matter. Because of its effect on structure formation, light 
scalar dark matter in the 10 -24 - 10 _2 ° eV range has been 
proposed [HI up to solve several long-standing astro- 
physical puzzles, such as the core-cusp, missing satellite, 
and too-big-to-fail problems [22] - DM self-interactions 
such as cubic or quartic potential terms likewise produce 
pressure contributions [33]; if they are sufficiently small 
in the early universe, they can be neglected in the present 
era J l3j . 

In Ref. [13!], it was pointed out that the interaction 
in Eq. 0 leads to a fractional oscillation in the fine- 
structure constant if <f) comprises the dark matter: 

a{t) ~ a [1 + d e K(j >o cos(w0f + 4)], (2) 

with a = e 2 /47r and 5 an arbitrary phase. The amplitude 
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(f>o depends on the ambient dark matter energy density 
Pdm ~ 0.3 GeV cm -3 as: 
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The amplitude <j> 0 scales as the square root of the en¬ 
ergy density, so even if (j) makes up a tiny fraction of 
the DM energy density, the oscillation in Eq. © may 
still be detectable. The field oscillation has an angular 
frequency equal to 77i0, and a coherence time of order 
2n/m ( j } v 2 , with v ~ 10” 3 the velocity dispersion of the 
DM in our Galaxy [23]. The amplitude in Eq. © scales 
inversely with the frequency of the oscillation, motivat¬ 
ing measurements of this effect at low frequencies. Our 
measurements are sensitive to angular frequencies below 
1 rad s” 1 , or DM masses 7710 < 10” 15 eV, for which 
we can take the DM oscillations to be coherent on the 
time scales of our experiment. At angular frequencies 
higher than 1 rad s” 1 , experiments with low noise floors 
such as gravitational-wave observatories and other reso¬ 
nant detectors may also reach sensitivity beyond existing 
bounds m, despite the smaller oscillation amplitude. 

Transition energies in atoms are convenient observables 
to look for low-frequency signals of changes in masses and 
couplings, and have been used to look for linear drifts in 
a [251 - 131] . Atoms are exactly reproducible systems across 
long time scales, their electronic structure solely deter¬ 
mined by masses and couplings of the Standard Model. 
While experimental imperfections can lead to fluctua¬ 
tions in observed transition frequencies, recent advances 
in optical metrology techniques have allowed for the de¬ 
termination of transition frequencies with stabilities and 
accuracies approaching 5vjv ~ 10” 18 [32H35] . 

In this work, we perform a spectroscopic analysis in 
two isotopes of dysprosium to search for the coupling in 
Eq. Q through the effect of Eq. 0. In the remainder 
of this Letter, we describe the atomic level structure of 
dysprosium and the experimental setup, followed by an 
account of the data analysis. No statistically significant 
oscillation consistent with a variation is observed; new 
limits are placed on couplings of DAI lighter than 10” 15 
eV. 

Dysprosium (Dy) is a rare-earth element with nuclear 
charge Z = 66 and 7 stable isotopes, of which we use 
those with atomic masses A = 162 & 164. A large num¬ 
ber of valence electrons endows Dy with a complex en¬ 
ergy level structure, including a nearly degenerate pair of 
opposite-parity states, denoted A and B for the even and 
odd states, respectively [37]. Spectroscopy of the radio¬ 
frequency (rf) electric-dipole transition between A and 
B revealed that their energy splitting corresponds to fre¬ 
quencies less than 2000 MHz, and that in some isotopes, 
A is the more energetic state ( 162 Dy), while in others 
( 164 Dy), B has the higher energy [33]. Figure [T] depicts 
the relevant energy level structure, with a focus on lb2 Dy 
and 164 Dy. 



FIG. 1. Energy level diagram of dysprosium showing the 
nearly-degenerate pair of states A and B, and the transitions 
used for state preparation and fluorescence detection. The 
inset shows the A and B states for 164 Dy (754 MHz) and 
162 Dy (235 MHz), and their response to a positive change in 
a: the transition frequency increases in 162 Dy and decreases 
in 164 Dy. The bottom-right graph displays a lineshape deriva¬ 
tive for the 754 MHz transition in 164 Dy. 


The energy splitting between the nearly degenerate 
pair is extremely sensitive to variation of the fine- 
structure constant. A change in a yields a frequency shift 
5v = ±v$ a 5a/a with v$ a « 2 • 10 15 Hz [33]. By compari¬ 
son, optical frequency measurements with trapped atoms 
and ions have similar absolute sensitivities to variation in 
a, but the near-degeneracy of the A and B levels relaxes 
requirements on the fractional accuracy and stability of 
the frequency reference. Spectroscopy in Dy resulted in 
one of the most stringent constraints on a present-day, 
linear variation of a at the level of |Aa/a| < 10” 16 yr” 1 , 
the best limit on a possible coupling of a to gravitational 
potential [25j , and stringent limits on violation of Lorentz 
invariance for electrons [3D] , 

A beam of dysprosium atoms is generated by an oven 
operating at 1400 K, and prepared in state B by two laser 
excitations followed by a spontaneous decay, as shown 
in Fig. [l] The lifetime of state B is long enough that 
it can be considered metastable in this experiment 331 . 
An rf electric held, whose frequency is compared with 
a cesium reference and a GPS-stabilized rubidium oscil¬ 
lator, then excites atoms to state A, which decays to 
the ground state through several channels, one of them 
emitting a 564-nrn photon. This fluorescence is detected 
with a photomultiplier tube (PMT), allowing the reso¬ 
nant frequency to be determined by maximizing huo- 
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FIG. 2. Simplified apparatus layout, a) Thermal beam of 
dysprosium atoms an oven heated to 1400 K. b) Skimmers 
collimate the beam, c) Linear polarizers suppress residual 
ellipticipy due to birefringence of vacuum chamber windows, 
d) Cylindrical lenses match the divergence of the atom and 
laser beams, for efficient population transfer BE- e) Mag¬ 
netic field coils compensate stray magnetic fields, f) Radio¬ 
frequency field interaction region, g) Photomultiplier tube 
for fluorescence detection. 

rescence with respect to the electric-field frequency. A 
simplified illustration of our setup is shown in Fig. [2] 
and more details can be found in Refs. [42] 03]. The 
optimal statistical precision is <j v = 7/(27 t\//Vt), where 
7/(27 t) ss 20 kHz is the natural linewidth of the transition 
(see bottom panel of Fig. [Tj) , N is the number of fluores¬ 
cent photons per unit time, and r is the total integration 
time. An estimate of N = 10 9 s -1 gives a statistical mea¬ 
surement precision of <j v \Jt = 0.6 Hz/-\/Hz. In the actual 
experiment, the statistical background is dominated by 
leakage into the PMT of blackbody radiation and scat¬ 
tered light from the 669-nm excitation laser, limiting the 
statistical precision for 162 Dy and 164 Dy to 4 Hz/-\/Hz. 
This corresponds to a sensitivity of 5a/a ~ 2 • 10 -15 
after one second of integration, and improves as 1 /yfr. 

Our data consist of transition frequency measurements 
on 162 Dy and 164 Dy—denoted by 1762 ~ 235 MHz and 
z/i 64 « 754 MHz, respectively—recorded at different 
times ti . The energy hierarchy of states A and B is 
reversed between these isotopes, which means that the 
response of their transition frequencies to a variation of 
a is opposite in sign, shown diagramatically in Fig. |Tj 
An increase in a would increase the transition frequency 
between the A and B states in 162 Dy, but decrease it in 
164 Dy. With measurements in both isotopes, two types 
of data sets can be constructed: one “out-of-phase” sig¬ 
nal data set maximally sensitive to a variation, the other 


an “in-phase” control channel data combination which 
can be checked for systematic effects. Because each iso¬ 
tope’s transition frequency was not observed simultane¬ 
ously, any systematic effects in the control data need to 
be subtracted to account for potential spectral leakage of 
these effects into the signal data. 

Measurements were performed on 10 different days 
over a two-year period (2010-2012), as in Ref. [25]. We 
split up the data into two different sets, one long-term 
data set (hereafter “LT”) of combined measurements on 
the first 9 measurement days, and a short-term data set 
(“ST”) of the last day of measurements. In the LT data 
set, the measurements of 7762 an d (764 are averaged over 
the course of each of the 9 days, such that each day- 
averaged measurement is dominated by a systematic er¬ 
ror of 0.48 Hz (1 Hz on the first day, due to a different 
experimental configuration). After subtracting the over¬ 
all mean in both LT sets, we define the LT signal data 
set as Signal {ti) = ^162 (^z) — ( 764 (^)j sensitive to vari¬ 
ations of a. Adding the frequencies from each isotope 
in phase defines the LT control data set ('control (ti) = 
(762 (ij) + ^i64(^i)- The ST data set comprises 2303 mea¬ 
surements of vis 2 and (764 taken over a span of 14.5 h 
on 19 October 2012. Data from each isotope are mean- 
subtracted, and combined into the ST signal time series 
data as Signal ({U,tj}) = (^({L}) - (764 ({tj})- Simi¬ 
lar to the LT data combination, we also construct a ST 
control data set ('control ({I7 tj})- 

Our aim is to search for harmonic variations in the 
dysprosium transition frequency data. A signal or limit 
at an angular frequency to can then be converted into a 
signal or limit for scalar dark matter with coupling d e 
and mass ~ w. We perform linear least-squares anal¬ 
ysis (LLA) on ('control (ti) and ^signal (ti) with waveforms 
of type v = Uq cos (cot + ip) + u c and v c a constant off¬ 
set m- The highest analyzed angular frequency w max 
is taken to be 47r/Af m i n , with f m ; n the shortest time be¬ 
tween measurements. LLA can also capture variations at 
angular frequencies smaller than the inverse time span of 
the data set T, when the waveform becomes essentially 
a frequency drift (near a node) or a quadratic frequency 
change (near an antinode). 

From the best-fit LLA waveform, we construct a nor¬ 
malized power spectrum: 

p (") - w 

where er 2 is the average variance of ('(ti), and Nq is 
the number of data points. For high angular frequen¬ 
cies 2n/T < oj < w max , the LLA power spectrum re¬ 
duces to the modified periodogram of Ref. U5|, and has 
simple statistical properties for white noise with equal 
uncertainties on each data point (a good first-order ap¬ 
proximation for both of our data sets). In this limit, 
the power P(w) at any frequency has a simple cumula¬ 
tive distribution function (CDF) which is w-independent: 
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FIG. 3. Power spectra of the best-fit waveforms for three com¬ 
binations of the ST data. The thin line and the dark and light 
bands correspond to the expected median and 68%- & 95%- 
CL intervals for the best fit. The raw control data set shows 
excess power at angular frequency cj s=s 1.4- 10 -3 rad s _1 (in¬ 
dicated by the arrow), which leaks into the raw signal data 
in nearby frequency bands. After fitting out this effect, the 
power in both data combinations is consistent with white 
noise, shown for the signal data in the bottom panel. The 
95% CL upper limit is depicted by a thick, dark blue line 
above the corrected best-fit signal power. 


Prob{P(w) < P} = 1 — e~ p . It follows that the expected 
median, and the 68% and 95% confidence intervals are 
0.69, [0.17,1.8], and [0.023,3.8], respectively. At inter¬ 
mediate and low angular frequencies uj < r/T, the sta¬ 
tistical behavior of the power is quite complicated, and 
we simulated its expected distribution with Monte Carlo 
(MC) using the same errors as on the data. 

We set a 95% CL upper bound on the power by as¬ 
suming the observed best-fit oscillation is a true signal, 
adding white noise via MC many times, and taking the 
95th-percentile value of the combined power. At high 
frequencies, an approximate, analytic expression for this 
value can be obtained (J5]. In the presence of a real 
signal with power P s , the CDF of the periodogram is 
modified to Prob(P(u;) < P \ P s ) = 1 — e~^ p+Pa ^(P, P s ), 
where 4>(P,P S ) = j2m=o XX=o P k P™/■ Given an 
observed power P s , we define our 95%-CL upper limit 
Piim by 0.95 « Prob(P(w) < Pii m |-P s )- At very low 
angular frequencies w <C 27r/T, where the periodogram 
statistics are invalid, the expected limit on the power 
(amplitude) scales as w -4 (w~ 2 ) because one cannot ex¬ 
clude being near an antinode of an oscillation. 

In Fig. [3] we plot the results of this analysis for the ST 
data set. Excess power is observed in the control data— 
sensitive only to in-phase fluctuations of vis 2 and v i 64 — 
at high significance: P(1.4-10~ 3 rad s _1 ) ss 57. We could 
not conclusively identify the source of this monochro¬ 
matic variation, observed independently in both isotopes 
with approximately equal amplitude, phase, and period 


FIG. 4. Upper limits at 95% CL on the coupling |d e | to pho¬ 
tons as a function of dark matter mass m^,. Shaded bands 
correspond to 68%- and 95%-CL intervals around the me¬ 
dian expected limit (thin line); thick lines depict the observed 
upper limits. Parameter space excluded at 95% CL by EP- 
violating force tests is shown in gray. 


of 4.3 • 10 3 s. Its best-fit amplitude of u 0 ss 0.78 Hz is 
comparable to the magnitude of a systematic error due 
to unstable electronic offsets [23], while its period may 
point to room temperature fluctuations. Variations due 
to reference clock instabilities are excluded because they 
would produce a larger effect in 164 Dy than in 162 Dy. 
Upon subtracting the best-fit waveform to the variation 
in the control data, which partially leaked into the signal 
data in neighboring frequency bands, the ST corrected 
power spectra were compatible with white-noise fluctua¬ 
tions. For the raw LT data, the observed power spectra 
in both the control and signal data were also consistent 
with white noise. 

Assuming the field 4> comprises all of the DM, power 
values can be translated to absolute values of the scalar 
coupling to photons via: 
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where ug a ss 2 • 10 15 Hz is the aforementioned sensitivity 
coefficient to fractional variations of a. Note that any 
bound on |d e | scales only as the square root of the energy 
density in the (f) field via Eq. ([3]). Figure [4] shows the 95%- 
CL upper limit on |d e | as a function of mass ~ w for 
both data sets. In particular, the observed limit from 
the ST data is the limit curve of Fig. [^transformed with 
Eq. ([ 5 ]). For masses m<f, < 3 • 10“ 18 eV, our ST limits on 
|d e | exceed those set by EP tests. In our most sensitive 
mass window—~ 1CU 22 eV, which coincides with the 
astrophysically motivated mass range |21l — we exclude 
DM couplings down to d e ~ 4.2 • 10~ 8 times gravitational 
strength. 

We have presented the best limit on neutral scalar dark 
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matter coupling to photons for masses below 3-10 -18 eV. 
It is also the first limit of its kind—exploiting the re¬ 
sponse of atomic transition energies to tiny fractional 
oscillations of the fine-structure constant—and we hope 
many similar searches will follow. Optical clock systems 
will likely improve on our limit of \d e \ with sufficient 
data. Combined with existing microwave and future nu¬ 
clear clocks, one can also achieve sensitivity to fractional 
variations in the electron, quark, and proton masses, ex¬ 
panding the potential discovery reach of light scalar dark 
matter to other couplings m- Small-scale precision ex¬ 
periments can thus search for cosmic dark matter with 
feeble interactions generated at the highest energy scales. 
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